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Abstract 

We consider the polynuclear growth (PNG) model in 1+1 dimen- 
sion with flat initial condition and no extra constraints. The joint 
distributions of surface height at finitely many points at a fixed time 
moment are given as marginals of a signed determinantal point pro- 
cess. The long time scaling limit of the surface height is shown to 
coincide with the Airyi process. This result holds more generally for 
the observation points located along any space-like path in the space- 
time plane. We also obtain the corresponding results for the discrete 
time TASEP (totally asymmetric simple exclusion process) with par- 
allel update. 



1 Introduction 

The main focus of this work is a stochastic growth model in 1 + 1 dimensions, 
called the polynuclear growth (PNG) model. It belongs to the KPZ (Kardar- 
Parisi-Zhang [19]) universality class and it can be described as follows (see 
Figure [T]). At time t, the surface is described by an integer- valued height 
function x t— > h{x,t) G Z, x G M, t G M+. It thus consists of up-steps (J") 
and down-steps (~L). The dynamics has a deterministic and a stochastic part: 

(a) up- (down-) steps move to the left (right) with unit speed and disappear 
upon colliding. 
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Figure 1: Illustration of the PNG height and its dynamics. The bold vertical 
piece is a nucleation. The arrows indicate the movements of the steps. A 
Java animation of the PNG dynamics is available at [9]. 



(b) pairs of up- and down- steps (nucleations) are randomly added on the 
surface with some given intensity. 

The up- and down-steps of the nucleations then spread out with unit speed 
according to (a). The PNG model can be interpreted in several different 
ways, see [11] for a review. 

On a macroscopic scale the surface of the PNG model grows determin- 
istically, i.e., lim^^oo ^~^^(^^, ^) = H{^) is a non-random function. However, 
on a mesoscopic scale fluctuations grow in time. This is called roughening 
in statistical physics and extensive numerical studies have been made [3]. 
Since the PNG model is in the KPZ universality class, the fluctuation of the 
surface height is expected to live on a t^^^ scale and non-trivial correlations 
are to be seen on a t"^^^ scale. Therefore, to have an interesting large time 
hmit, we have to rescale the surface height as 

One of the most natural initial conditions for PNG is the fiat initial 
condition, i.e., h{x,Q) = for all x G M. We consider nucleations occurring 
with translation-invariant intensity^ We refer to the PNG model with such 
initial condition as flat PNG. In this case, by mapping the fiat PNG to a 
point-to-line directed percolation model it was proven [2,21,22] that the one- 
point distribution is, in the t — > cxo limit, the GOE Tracy- Widom distribution 
Fi, first discovered in random matrix theory [27]. However, no information 
on joint height distributions at more than one point has been previously 
known. 



""^In other words, the nucleation events form a Poisson process with constant intensity 
in the space-time upper half-plane. 
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New Results. The main results of this paper are precisely the computation 
and asymptotic analysis of these joint distributions. In particular, we prove 
the convergence of the height rescaled as in (11.11) to the Airyi process in the 
t ^ oc limit (see Section 12.21 for a definition of the process). The Airyi 
process has been discovered in the context of the asymmetric exclusion pro- 
cess [5-7,24]. Our result, stated in Theorem El is obtained by first determin- 
ing an expression for the joint distributions for finite time t (Proposition H]) 
and then taking the appropriate scaling limit. 

Proposition H] is actually just a particular case of Theorem [5l where we 
determine joint distributions along any space-like paths (as in Minkowski 
diagram), for which fixed time is a special case. The scaling limit is analyzed 
at this level of generality, thus Theorem [U] holds for any space-like paths. In 
contrast to previous works on the subject, our approach does not rely on the 
so-called RSK correspondence (RSK for Robinson-Schensted-Knuth), which 
was successfully applied for corner growth models, but does not seem to be 
well suited for the flat growth. 

On the way of getting the results for the flat PNG, we consider a dis- 
crete time version of it, the Gates- Westcott dynamics [12,23]. This model is 
closely related to the totally asymmetric simple exclusion process (TASEP) 
in discrete time with parallel update and alternating initial conditions. The 
corresponding results for this model are Theorem [T] for the joint distributions 
along space-like paths, and Theorem [3] for the convergence to the Airyi pro- 
cess in the scaling limit. For the TASEP, the extreme situations of space-like 
paths are positions of different particles at a fixed time and positions of a 
fixed particle at different time moments (tagged particle). The space-like 
extension for TASEP is based on the previous paper [4]. 

Previous works on PNG. Another type of initial conditions for the PNG 
model has been analyzed before. It is the corner growth geometry, where 
nucleations occur only inside the cone < t}. The limit shape if is a 
semi-circle, and the model is called PNG droplet. In this geometry, the limit 
process has been obtained in [23]; it is known as Airy2 process (previously 
called simply Airy process). The approach uses an extension to a multilayer 
model (inherited from the RSK construction), see [16,23]. The multilayer 
method was also used in other related models [8,13,14,17,18,26]. Also, 
for the fiat PNG it was used to connect the associated point process at a 
single position and the point process of GOE eigenvalues [10]. Results on the 
behavior for the PNG droplet along space-like paths can be found in [8] . For a 
very brief description of the previously known results on TASEP fluctuations 
see the introductions of [4, 25] . 
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Outline 



In section |2], we introduce our models and state the results. In section [HI we 
give an expression of the transition probability of the discrete TASEP as a 
marginal of a determinantal signed point process. In section H] the Fredholm 
determinant expression for the joint distributions is obtained. The argument 
substantially relies on the algebraic techniques of [4]. In sectional we consider 
the scaling limit of the parallel TASEP. In section [HI the continuous time 
PNG model is considered. In section U\ we consider the scaling limit for the 
continous PNG model. 
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2 Models and results 

We start from the discrete time TASEP with parallel update. Then we will 
make the connection with a discrete version of the PNG, from which the 
continuous time PNG is obtained. 

2.1 Discrete time TASEP with parallel update 

We consider discrete time TASEP with parallel update and alternating initial 
conditions, i.e., particle i has initial position Xi{0) = —2i, i e Z. At each time 
step, each particle hops to its right neighbor site with probability p = 1 — q 
provided that the site is empty. The particle positions at time t is denoted 

by Xi{t),i G Z. 

The dynamics of a particle depends only on particles on its right. This 
fact allows us to determine the joint distributions of particle positions also for 
different times, but restricted to " space-like paths" . To define what we mean 
with "space-like paths", we consider a sequence of couples (ni,tj), where rii 
is the number of the particle and ti is the time when this particle is observed. 
On such couples we define a partial order -<, given by 

{rii.ti) -< {nj,tj) if Uj > ni,tj < ti, and the two couples are not identical. 

(2.1) 
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A space-like path is a sequence of ordered couples, namely, 

<S = {{nk,tk),k = 1,2, . . . \{nk,tk) -< (rifc+i, tfc+i)}. (2.2) 

The reason of the name "space-like" will be clear in the large time limit, 
where everything becomes continuous. Then space-like is the same concept 
as in the Minkowski diagram. The border cases for space-like paths are fixed 
time {ti = t,\/i) and fixed particle number (n^ = n,\fi). The next theorem 
concerns the joint distributions of particle positions. 

Theorem 1. Let particle with label i start at Xi{0) = —2i, i G Z. Consider 
a space-like path S. For any given m, the joint distribution of the positions 
of the first m points in S is given by 

m 

P( n {^"ki^k) > flfc}) = ^^^(1 - X^^^xi"^)f2({(„i,ti),...,(„„,t„)}xz) (2.3) 

k=l 

where X^c^\{^k^tk)i'^) = 1(3^ < ct/c)- The kernel Kt is given by 

Kiim,t,),X,; in2,t2),X2) = -0«"l'*l)'("^'*^))(a^l, a;2)l[(niA)^(n„t.)] 

+ ^((ni,ti),xi; (712,^2), xs), (2.4) 

where 

K{{ni,ti),xi; {n2,t2),X2) 

-1 /• ^^(1 + z)^'2+'^i+"2 ^i_pyi~2m-xi 

27ri Jp ^__2^a;i+ni+7i2+l _|_ p2^ti+t2+l-(a;i+ni+?i2) ' V • / 

and 

1 / , (l+pwY^'^^ 

= — (h aw- r — 

27ri /p_^ (1 + w)^'i-^'2+i 

where Tq (resp. is any simple loop, anticlockwise oriented, with (resp. 
—1) being the unique pole of the integrand inside the contour. 

Remark. In the limit p — > under the time scaling by p~^ the discrete 
time TASEP converges to the continuous time TASEP, and Theorem [T] turns 
into a special case of Proposition 3.6 of [4], where a more general continuous 
time model called PushASEP was considered. 



—w 



(1 + w){\ + pw) 



"1 — "2 



(2.6) 
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2.2 Airyi process and scaling limit 

Starting from Theorem [T] we can analyze large time limits. The limit process 
is the so-called Airyi process introduced in [6,24], which we recall here. 

Definition 2 (The Airyi process). Define the extended kernel, 

^.Ai(ri,6;^2,6) = / exp (-^ — ^) 1{t2 > Ti) 

v/47r(r2 - ri) V 4:{t2 - n) J 

+Ai(6 + 6 + (^2 - nf) exp (^(r2 - ri)(6 + 6) + ^(^2 - ri)^^ . (2.7) 

r/ie Airyi process, Ai, is the process with m-point joint distributions at 
Ti < T2 < . . . < given by the Fredholm determinant, 



p(r\{Mn) < Sk}] = det(l - XsKa, 



Xs)l2({^i,...,^„}xr), (2.^ 



where Xs{n,x) = l{x > Sk)- 



Theorem [T] allows us to analyze joint distributions of particle positions 
for situations spanning between fixed time and fixed particle number (the 
tagged particle problem). One way to parametrize such situations is via a 
space-like path. We thus consider an arbitrary smooth function vr satisfying 

|7r'(^)| < 1 and 7r(^) + > 0, (2.9) 

see Figure [21 The requirement 7r{6) + > reflects t > 0. Then, we 
choose couples of (t, n) on {((7r(^) + e)T, {it (6) - e)T), 6 eR}, where T is a 
large parameter. The case of fixed time, say t = T, is obtained by setting 
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n{9) = 1 — 9, while fixed particle number, say n = aT, by tt{6) = a + 9 with 
some constant a. 

From KPZ scaling exponents [19], we expect to see a nontrivial limit if 
we consider positions at distance of order T^/^. Thus, the focus on the region 
around 9T is given by 9T — uT'^^^, i.e., setting 9 — uT~^/^ instead of 9 and, 
by series expansions, we scale time and particle number as 

t{u) = [{n{9) + 9)T - {tx\9) + l)uT''l^ + \tx"{9)u^T^'^\ , 
n{u) = [(7r(e) -e)T + (1 -7r'(^))MT2/3 + l7r"(^)u2Ti/3j. (2.10) 

The KPZ fiuctuation exponent is 1/3, thus we expect to see fiuctuations of 
particle positions on a scale of order T^/^. Therefore, we define the rescaled 
process by 

Xn(^)(t{u)) - {-2n{u) + wt{u)) 
u^^t[u) = ^^ -j^^ . (2.11) 

Here the mean speed of particles, v, is determined to be v = 1 — from 
the subsequent asymptotic analysis but can be known beforehand from the 
stationary measure for density 1/2 [15]. This process has a limit as T — oo 
given in terms of the Airyi process. 

Theorem 3. Let he the rescaled process as in 1^2. Then 

lim Et{u) = K^Ai^Khu), (2.12) 

T — ^oo 

in the sense of finite dimensional distributions. The vertical (fluctuations) 
and horizontal (correlations) scaling coefficients are given by 

= {tx{9) + 9f'\l - qfl\^/\ (2.13) 
(7r(g) + g)^/3(l-g)^/3g-i/6 
(7r'(e) + l)(l-Vg)/2 + l-7r'(^^)- ^ " ^ 

Remark. A similar result for the PushASEP with alternating initial 
condition has been proved in Theorem 2.2 of [4]. 



2.3 TASEP and growth models 

As mentioned in the introduction, the discrete TASEP with parallel update is 
related to a surface growth model from which the polynuclear growth model 
in continuous time can be obtained as a limit. Let t > and x e M denote 
the time and the one- dimensional space coordinate respectively, and let ht{x) 
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Figure 3: A surface growth model. For half-odd integer times this is equiva- 
lent to the discretized Gates- Westcott dynamics and for integer times to the 
discrete TASEP. 



be the height of the surface at time t and at position x. Let us introduce 
a dynamics of ht{x) as follows. Initially, at time t = 0, the surface is flat; 
ho{x) = 0, for all x G M. Right after each integer time (t = 0+, 1+, 2+, . . .), 
there could occur a nucleation with width and height 1 with probability q 
(0 < g < 1) independently at each integer position x such that t + x + ht{x) is 
even. Each nucleation is regarded as consisting of an upstep and a downstep 
and each upstep (resp. downstep) moves to the left (resp. right) with unit 
speed. This is a deterministic part of the evolution. When an upstep and a 
downstep collide, they merge together. See the solid line in Figure [3] for an 
example until t = 2. The dynamics of the growth model, if we focus only 
on half-odd times (t = |, |, . . .), is the same as one considered in [23], i.e., 
a discretized version of the Gates- Westcott dynamics [12]. It is known that 
in an appropriate g — limit this growth model reduces to the standard 
continuous time PNG model [23]. 

To see the connection to the discrete TASEP, let us focus on integer 
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times {t = 0, 1, 2, . . .) and positions {x G Z) from now on and represent the 
surface as consisting of elementary upward slopes / and downward slopes 
\ as indicated by dashed lines in Figure [3l At t = 0, even (resp. odd) x's 
are taken to be the center of the upward (resp. downward) slopes. Then 
the dynamics of the surface is described as follows: At each time step the 
surface grows upward by unit height deterministically and then each local 
maximum (/X) of slope turns into a local minimum (\/) independently 
with probability p = 1 — q. If we interpret an upward (resp. a downward) 
slope as a site occupied by a particle (resp. an empty site), this is equivalent 
to the discrete time TASEP with parallel update under the alternating initial 
condition. 

The relation between the surface height ht{x) and the position of the 
TASEP particle is given by 

htix) <H^ x^ t-.-H j it) > X (2.15) 

and is understood as follows. On the plot of the surface at some fixed time 
t, draw also the initial surface aX h = t. See Figure H] for an expample. 
Then, from the correspondence between the growth model and the TASEP, 
the surface at time t can be regarded as the particle positions. In this plot 
particles move along the down-right direction as indicated. The left hand side 
of (12.151) is equivalent to the condition that the TASEP particle corresponding 
to (x, h = t~H) has already reached x. Since the axis of the particle number 
n is in the down-left direction, the value of n corresponding to {x,h = t — H) 
is \ {t — H — x)/2\. This consideration results in the relation (I2.15p . From 
the relation (I2.15P the joint distributions of the height of the growth model 
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is readily obtained through 

m m 

p[[]{hu{x^)<H,}^=p(^f]{xnXU)>x^}), (2.16) 

i=l i=l 

combined with Theorem [H 

When q ^ 0, the TASEP particles move almost deterministically and 
the surface ht{x) grows slowly, when a particle decides not to jump (with 
probability q). The continuous time PNG model is obtained by taking g — > 
while setting space and time units to ^/q/2 (the 2 is chosen to have nucleations 
with intensity 2 like in [23]). Denote by x and t the position and time 
variables in the continuous time PNG model. The PNG height function 
/i™^(x, t) is then obtained by the limit 

/i™«(x,t) = lim/i2t/v^(-2x/V^). (2.17) 

Here the minus sign on the right hand side is put for a convenience. The 
results below do not depend on this sign because of the symmetry of the 
model in consideration. The joint distribution of the surface height at time 
t is given as follows. 

Proposition 4. Consider m space positions xi < X2 < . . . < x^- Then, the 
joint distribution at time t of the heights /i^^^(xfc, t), k = 1, . . . ,m, is given 
by 

m 

P( fl {/i™«(xfc,t) < Hk}^ = det(l - XHKr''xH)iH{.„...,.^}xZ) (2.18) 

k=l 

where the kernel is given by 

ir™^(Xi, hi- X2, h^) = -I\h,-h,\ (2(X2 - Xi)) 1(X2 > Xi) 

+ ( 2t-X2 + x! j '^'^^"^ (2v/4t2-(x2-xi)2) l(2t > |X2 - xi|) 

(2.19) 

where In{x) and Jn{x) are the modified Bessel functions and the Bessel func- 
tions, see e.g. [1]. 

The last indicator function is obvious if one thinks about the PNG model. 
In fact, the height at position x at time t depends on events lying in the back- 
ward light cone of (x, t) on M x IR+. Thus, when |x2 — xi | > 2t, the backwards 
light cones of (xi,t) and (x2,t) do not intersect in M x R+, which implies 
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that the two height functions are independent. The Fredholm determinant 
then sphts into blocks. 

The result of Proposition H] is actually a specialization of a more general 
situation which follows from the TASEP correspondence. In the TASEP, the 
space-like paths vr we had for particle numbers and times become the paths 

(x,t) = (7r(0) -30,0 + 7r(^)). (2.20) 

The condition |7r'(^^)| < 1 implies that 9t/9x G [—1,0], i.e., these are space- 
like paths as in special relativity oriented into the past. By the symmetry of 
the problem, one can consider also space-like paths locally oriented into the 
future, just looking at the process in the other direction. 

Denote by 7 such a path on M x M.*^, i.e., (x, t = 7(x)), then 6 and 7r(0) 
are given by the relations 

9 = (7(x) - x)/4, n{9) = (37(x) + x)/4, (2.21) 

and the joint distributions of the surface height along the path 7 are expressed 
as in Theorem [5l 

Theorem 5. Let us denote by tk = 7(xa:). Then, the joint distributions of 
/i^^^(xfc, tic) , k = 1, . . . ,m, is given by 

m 

P( fl {/i™«(xfc,t,) < Hk}^ = det(l - XHi^™''x//)^2{{(x,tO,...,{x.,t,)}xZ) 

(2.22) 



k=l 



where the kernel is given by 
ir™«((x„tO,/.i;(x„t,),/.^- I""' ^^ + '^ 



X2 — Xi — ti + t2 ^ 

X ^\hi-h2\ (2 a/ (x2 - Xi)2 - (t2 - ti)2j l{(ti+xi,ti)^(t2+X2,t2)} 

/ (ti+t2) + (x2-x0 \ '^'^^'^^Z' ( X 

+ Ut, + t2)-(x2-x0 J -^'^^'-^ V^^^' + - - j 

X l(ti +t2 > |xi -X2I) (2.23) 

where In{x) and Jn{x) are the modified Bessel functions and the Bessel func- 
tions. The condition I{(ti+xi,ti)^(t2+x2,t2)} means that X2 — xi > ti — t2 > 
or X2 — Xi > ti — t2 > 0. 

In the first term, for X2 > Xi, the condition X2 — Xi > ti — t2 is satisfied 
for tk = 7(xfc). Also, notice that when X2 — xi ^ ti — t2, the first term of 
the kernel goes to (2(x2 - xi))l'^i-'^2|/(|/^^ _ /^^l)! . 
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Figure 5: A space-time path 7 for continuous time PNG. T is proportional 
to the PNG time t. 



2.4 Scaling limit for the continuous PNG model 

The last result of this paper is the large time behavior of the flat PNG. The 
large parameter denoted by T is proportional to time t. Using the function 
7, we consider t = T7(x/T), see Figure O 

Since the system is translation invariant, we focus around the origin, i.e., 
we look at the PNG height at 

t{u) = 7(o)r + y(o)MT2/3 + 17"(o)m2ti/3. ^ ■ ^ 

The surface height grows with the speed equal to 2. Thus, for large time t, 
the macroscopic height will be close to 2t. Fluctuations live on a T^/^ scale. 
Consequently, we define the rescaled height process by 

(2.25) 

The large T (thus large time too) behavior of /i™^ is given in terms of the 
Airyi process as stated below. 

Theorem 6. Let /i™^ be the rescaled process as in ^2.25\) . Then, in the 
limit of large T, we have 

lim h™^{u) = S^AiiShu), (2.26) 

T — >oo 

in the sense of finite dimensional distributions. The scaling coefficients 
and Sh are given by 

= (27(0))l/^ Sh = (27(0))^/^ = Sl (2.27) 
For 7(0;) = 1, i.e., fixed time, this was conjectured to hold in [6]. 
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3 Transition probability for the finite system 



Let G{xi, . . . ,xi\i;t) denote the transition probability of the parallel TASEP 
with particles starting at t = at positions < ■ ■ ■ < yi- This is the 
probability that the particles starting from positions y^ < . . . < yi at 
t = are at positions xn < ■ ■ ■ < Xi at t. 

Consider a determinantal signed point process on the set 
X = {x^, 1 < i < n < N} by setting the measure 

Wn{x) = I Yl det(0«(x^,x^+|))o<ij<„ j det{F_i+i{xf-yN+i^i,t+l-i))i<ij<N 

^ n=l ' 

(3.1) 

where 

{1, y>x, 
p, y = x-l, (3.2) 
0, y<x-2, 

the function F„(x, t) defined by 

If -jy" 

FJx, t) = (p dw- ^^(1 + (1 - q)wy, (3.3) 

^ ' ^ 27ri/ro_i (l + w)"+^+i^ ^ ^ ' ^ ^ 

and where we used the convention, Xq = — oo. 

The following proposition states that the one time transition probability 
of the TASEP is given as a marginal of the signed measure (13.11) . 

Proposition 7. Let us set = Xn, n = 1, . . . ,N. Then 

G{xi,...,XN;t) = Y,WN{x) (3.4) 

V 

where summation is over the variables in the set, 

V = {x^, 2<i<n< Nlx'l > x[Li} (3.5) 

varying over Z. 

Note that W]y{x) is actually symmetric with respect to permutations of 
variables with same upper index, so the ordering in (I3.5jl is used for singling 
out the minimal x^ = min{x", i = 1, . . . ,n}. 

Remark. Similar representations for the transition probability of con- 
tinuous time TASEP, discrete time TASEP with sequential update and 
PushASEP have been obtained in [4-6]. 

In the different parts of the proof of Proposition U\ we will use several 
properties of the function F„, which are listed below. 
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Lemma 8. 

oo 

F^+,{x,t) = Y,F'n{y,t), (3.6) 



y=x 



F„(x, t + l) = qF„{x, t) + (1 - g)F„(a; - 1, t) (3.7) 

= F„(x,t) + (l-g)F„_i(x-l,t), (3.8) 

(</.»*F„)(x,t) = F„+i(x,t + l), (3.9) 

F_„(x, -n) = for X < -n, n > 0, (3.10) 

Fn{x,n) = for X > n, n > 0, (3-11) 

F„(n,n) = (l-g)",r2>0, (3.12) 

F_„(-n, -n) = l/(-g)", n > 0. (3.13) 

iiTere represents the convolution: (0" * /)(2;) = <P\x,y) f {y) . 

Proof of LemmalB These are proven by using the definition (13.21) and (13.31) . 

□ 

The first step in the proof of Proposition [7] is the following Lemma. 
Lemma 9. Let us set 

{yV-^ ^ y > X, 

1-q, y = x-l, (3.14) 
0, y<x-2 

and 0f,(— oo, y) = . Then, for any antisymmetric function /(6i, . . . , bn), 
^ det(0^(ai, bj))o<i,j<n ■ /(&i, • • • , ^n) 

>...>fei>feo 
bf): fixed 

= guiai,bo) ^ det(0|(ai,6j)) 

l<i,j<n f{b,,...,bn) (3.15) 

>...>fei>feo 
bo: fixed 

where On > . . . > ai, = — oo and 

'O, b>a, 

g,ia,b) = {u\l-{l~q)u), b = a-l, (3.16) 

b<a-2. 
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Proof of Lemma From the antisymmetry of / and of the determinant, 
(I3.15P is equivalent to 



0<i,j<n n) 



bi,...,b„>bo 
bo: fixed 



f{b 



li ■ ■ ■ 1 J 



(3.17) 



6i,...,6„>6o 
6o: fixed 



Since a basis of the antisymmetric functions is made of the antisymmetric 
delta functions and the relation to prove is linear in /, it is enough to consider 



f{bi, ...,bn) 



{-ly, if (6i, ...,bn) = (&ai, • • • , b^J for some a e Sn 



0, 



otherwise 



(3.18) 

for fixed > bo. Here Sn is the group of all permutations of 

{l,...,n}. For this special choice of /, the left hand side of (13.171) is n! 
times the single determinant, 



det 



,bo 



(3.19) 



^lia-nM) 0*(a„,6i) ... 0|(a„,6„) 

We have the following three cases, 
(a) cti < &o- the second row gives {]/"^~°'^ ^ . . . , y^^-°-^'^ which is proportional 

to the first row. Therefore in this case the LHS is zero. 

(b) ai = 6o + 1- The second row is (1 — g, v^^~°-^ . . . , z/^"~"i). Subtracting 
v°-^ times the second row from the first row one obtains 



z/^«(l - (1 - q)v) ■ det(0«(a,,6,))i<,,,<„. 



(3.20) 



(c) ai > &o + 1: The first column is (z/''°, 0, . . . , 0)*. Thus the determinant is 
z/^o ■ det(0j,(ai,6j))i<ij<„. 

The result in each case agrees with n! times the RHS of (13.171) and hence the 
lemma is proved. □ 

Let M{xi, . . . ,xn) denote the number of j's s.t. xj — ^j+i — 1)J — 
1, . . . , — 1. Using the above lemma with z/ = 1 in which case 0j, reduces 
to (p^, we have the following result. 
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Lemma 10. With = x„, n = 1, . . . , N , one has 

Wn{x) = det[Fj_i{xN-j+i - yN-i+i,t + 3 - i)]i<r,j<N 

(3.21) 



V 



— : G{xi, . . . , XN',t). 
Proof of Lemma [73 For simplicity, we denote 

fi{x) = - yN-^+l, t-i + 1), (3.22) 

for i = 1, . . . , N. From the definitions 03.11] . the LHS of fl3.2ip writes 

J2 [ nd^*(<^na^r,a:Si'))o<M<n)det(/.(a;f))i<,,,<^. (3.23) 



<>a;;i_i>...>a;5^ ^ n=l 
x":fLxedA<n<N 



Applying Lemma [9] with u = l,n = N — 1, ai = x 



bi = xf^i, i = 0, . . . , N — 1 and 



Af-l 
i 1 



/(&1, ...,bn) = det{fi{Xj ))i<i,j<N, 



1 iV- 1, 



(3.24) 



we obtain 

m^=gi{x^-\x^] 



.N-2 . 



a;;i><_i>...>x!|' ^ n=l 
x":fLxedA<n<N-l 



X 



det(0''(xf \x^i))i<ij<Ar_i ■ det(/i(xf ))i<ij<Ar. 



X fixed 



(3.25) 
Heine's identity, 

^ ^ det((/)i(xj))i<ij<„det(V'j(xj))i<ij<„ = det [0i * ^j]^^,.^^.^^, (3.26) 



allows us to rewrite the last summation in f l3.25p as 

(0«*/i)(xf"^) ... (0«*/o(xj^:}) 

/^(a:f) (0«*/^)(xf-i) ... (0«*/^)(xj^:}) 



det 



(3.27) 
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We repeat the procedure up to a total of j — 1 times in column j and we get 



(!325D= ( n^?i(^^^^') det[0«*...*0«*/,(a;f-^+i)]i<,,,<^. (3.28) 

^ n=l ^ 

The proof of the lemma is finished using fl3.22p . (13.91) and 
Un=i 9ii^i^^i'^^) = □ 

Proof of Proposition We need to prove 

G{xi, . . . , xn; t) = G{xi, . . . , xn; t). (3.29) 

This statement was also proved in [20] by the Bethe ansatz techniques. Our 
proof is by induction in t. We start by showing that the initial conditions 
agree, i.e., G{xi, . . . , xn] 0) = G{xi, . . . ,xn;0), that is, 

N 

^(xl,...,x^) . (iet[Fj_i{xN-j+i - yN-i+i,J - i)]i<ij<N = n ^^n,yr.- (3-30) 

n=l 

We first show that LHS of (I3.30p is zero if xn Vn- If xn < Un — 1, since 
UN^i+i > Vn + i - one has xn - VN-i+i < -i + l,i = 1,. . .,N. Then, 
from (13.101) we have Fi_i{xN — Vn-i+i, ^ — i) = 0, i.e., the first column of 
LHS of (I3.30p is zero. Similarly, ii x^ > i/n + 1, since xtv-j+i > + j — 1, 
one has x^-j+i — Vn > j — = 1, . . . , A^. Then, from (13. lip we have 
Fj_i(x7v-j+i — yNii — 1) = 0, i.e., the first row of LHS of (I3.30p is zero. This 
agrees with RHS of (I3.30p also being zero if xn Vn- 
Now let us assume xn = Vn- There are two cases. 

(a) vn-i > l/AT + l. In this case, since XN-yN-i+i = VN-VN-i+i < + = 
2, . . . ,N, one has Fi_i{xN — i/N-i+i, ^ — i) = 0,i = 2, . . . , N. Then the first 
column of LHS of (I3.30p is (1, 0, . . . , 0)* and hence the determinant is equal 
to det[Fj_i{xN-j+i - VN-i+iJ - i)]2<i,j<N- 

(b) Un-i = + 1- First let us see that LHS of (I3.30p is zero when xjy^i ^ 
Vn-i- We have xn~i > + 1 = + 1 = yN-i- If xn~i > Vn-i + 1, we 
have XN~j+i - Vn > xn-i + j - 2 - {vn-i - 1) > j, for j = 2,..., N, and 
XN-j+i — Vn-i > j — 1, for j = 2, . . . , A^. Then the first and the second 
row of LHS of (I3.30p are both of the form, (*, 0, . . . , 0) where * represents an 
arbitrary number and hence the determinant is zero. Hence LHS of (I3.30p is 
zero if xjy^i ^ Vn-i- On the other hand, when x^-i = Vn-i, the upper-left 
2x2 submatrix of the determinant is 



Fo(0,0) Fi(l,l)" 




" 1 1-q 


F_i(-1,-1) Fo(0,0)_ 




-1/q 1 



(3.31) 
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whose determinant is 1/q. 

Repeating the same procedure, at each step one has either case (a) or (b). 
The final result is that = x^, for A; = 1, . . . , A^, otherwise the determinant 
in LHS of f l3.30p is zero. Moreover, when yk = Xk, k = 1, . . . , N, denote by 
rii, ni + n2, . . . ,ni + . . . + rii the values of j such that xj-i — Xj > 1. Then 
LHS of fl3.3Up is equal to 11™=! -^"m with 

D„ = det [Fj.iU -I, J- 0]i<,,,<, (3.32) 

Finally using (13.121) . (13.131) . we obtain an explicit form of the matrix. To 
compute its determinant it is enough to develop along the first row. The 
determinant of the (1, 1) minor is -Dn-i? while the one of the (1,2) minor is 
(— l/g)Z)„_i because the minor is the same as the (1,1) minor except the 
first column is multiplied by —l/q- AH the other minors have determinant 
zero, because the first two column are linearly dependent. Thus, Dn = 
1 ■ Dn-i — (1 — q)/{^Q)Dn-i, and since Di = 1, it follows that 

Dn = (3.33) 

This ends the part of the proof concerning initial conditions. 

Next we prove that (13.291) holds for t + 1 if it does for t. Since this is true 
for t = 0, by induction it will be true for all t E N. G satisfies the TASEP 
dynamics, thus 

G{xi, . . . ,XN]t + 1) 
= ^G{zi,...,ZN,t)w{z,x) = ^G{zi,...,ZN,t)w{z,x) (3.34) 



^w;(2;,x)g^("^--'"^Met[F,_i(2;^_,+i - + J - i)]i<i,. 



j<N- 



Here 



A'' I -I-) 1; Xn Z^i 

W{Z,X) = ^Vn, Vn=<^q, Zn < Zn-l - I, Xn = Zn, (3.35) 

^ 1; -^n ~l~ 1; 



n=l 



and in the second equality we have used the assumption of the induction. 
We rewrite (^(xi, . . . , xat; t + 1) using (13. 7p and (13. 8p as follows. For k from 
1 to N: 

(a) if Xk = Xk+i + 1, then we use (13. Sp to the + 1 — fcth column. Then, the 
new term with the (1 — Q') factor in front cancels out because it is proportional 
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t + 1 



t + 1 




q-1 



q-1 



(d) 



Figure 6: Graphical representation of (13.381) . The dots represents empty 
places, while a line leaving/ arriving to a point is an occupied position. In (a) 
and (b), on the left (resp. right) we indicate the weights different from 1 of 
LHS (resp. RHS) of (13.381) . In (c) and (d) the bottom and top lines of two 
blocks at distance 2 at time t+1 are represented, for the cases corresponding 
to (a) and (b) for the top block. 



to its left column of the determinant. 

(b) if Xk > Xk+i + 1, then we just use (13. 7p . 

With these replacements we get 

G{xi, . . . ,XN',t + 1) (3.36) 
= ^w{z,x)q^^^''---'''''^ det[Fj_,{zN-j+i - yN-i+i,t + J - i)]i<i^^^^ 

z 

where 

JV I 1, Xji 3Jn+l ~l~ 1; ■^n -^ni 

W{Z, X) =Y\ Vn, Vn= Iq, Xn > Xn+1 + 1, Z„ = X„, (3.37) 

^1 Q', Xji > Xji-i^l ~\~ 1, Zn Xji 1. 

Comparing (I3.34p and (I3.36p . it is enough to show 

g^("i'-'^^)w(z, x) = g^("i'-'"^)w(z, x). (3.38) 

This indeed holds and can be seen by checking case by case. We illustrate 
it using Figure El First consider a block of particles, say m of them at time 
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t + 1. There are two possibility of reaching this situations in one time step, 
as indicated in Figure E] (a) and (b). The products of all the weights on 
the right and on the left are the same, i.e., (13.381) holds for a single block 
of particles. If two blocks of particles at time t + 1 are at distance at least 
2, they are independent during one time step. We just have to check that 
fl3.38p holds for two blocks at distance 2 at time t + 1. Case (a) is illustrated 
in (c) and the weights are unchanged for both blocks. Case (b) is illustrated 
in (d). This time, the q on the top line of the second block becomes a 1, but 
this is compensated by an extra factor q on the left. □ 



4 Joint distributions along space-like paths 

Theorem 11. Let us consider particles starting from yi > y2 > . . . and 

denote Xj{t) the position of jth particle at time t. Consider a sequences of 
particles and times which are space-like, i.e., a sequence of m such couples 

5 = {{nk,tk), k = 1, . . . , m|(nfc, tfc) -< ('^fc+i, ^fc+i)}- The joint distribution of 
their positions Xn^{tk) is given by 

m 

P( n {^"fc(^fc) ^ "fe}) = - Xai^Xa)£2{{KA),...,{n™,im)}xZ) (4.1) 

k=l 

where Xa{{nk,tk),x) = l{x < a^). Here K is the extended kernel with entries 



n2 



k=l 



(4.2) 

where 

0*«"^'*i)'("2'*2»(xi,X2) (4.3) 
- 1,.. ^^(1 + ^)---+! [{l + w){l+pw)) 



and the functions "^^-i '^'^^ given by 



n—k 

■ (4.4) 



^-\ix) = ^i d^^i±^r 1 

^-^^ ' 27ri/r„_^ {l + wY-y^+^\{w + l){l+pw)) 

The functions are determined by the two conditions 

*::!,(^)C-fc(^) = ^KU l<k,l<n, (4.5) 
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and 



span{$":^^(x), k = 1, . . . ,n} = span{l, x, 



}• 



(4.6) 



The paths Fq,-! are any simple loops anticlockwise oriented including 0,-1 
and no other poles. 

Proof of Theoremllli This is the analogue of Proposition 3.1 of [4]. The first 
step to obtain this proposition, was Lemma 4.4 of [4]. In the simple case 
where in the model studied in [4] all particles have the same jump rates, i.e., 
Vi = 1 for all i, then Lemma 4.4 of [4] and our Proposition [7] have exactly the 
same structure. For the comparison of the two, we need just the following 
identifications. 



PushASEP 


Parallel TASEP 


Fn{x,a,b) 


(f)^{x,y) 

Fn{x,t) 



where 



Fn{x,t) = Fn{x,t + n). 
From 03.91) it then follows that F^ satisfies 

(0»*F„)(x,t) = F„+i(x,t). 



(4.7) 



(4.^ 



The identity (14.81) corresponds to Lemma 4.3 of [4]. To obtain Proposition 
3.1 of [4] one needs to prove Lemma 4.5 of [4] which needs as ingredients 
only Lemma 4.3 of [4]. These Lemmas are the only ingredients used to 
obtain Theorem 4.1 and then Theorem 4.2 of [4]. This last theorem concerns 
some general determinantal measures (a generalization of Lemma 3.4 in [6]), 
whose specialization to the PushASEP model is Proposition 3.1 of [4]. All 
these steps go through unchanged in our case provided we make the above 
identifications. 

To obtain the precise expressions of Theorem [TT], we use an integral rep- 
resentation for 0", namely, 

1 {l+pw){l + w) 



dw 



{w + l)^-f+i 



w 



From this we get 



1 

27ri 



dw- 



1 



Also, the Tt.^t-_^{x,y) from [4] now becomes 

1 

27ri 



'I +pw){l + w) 



w 



Fo{x -y,ti- ti 



i-i 



dW' 



[1 + pwf 



'1 + w)'-^-y+^ 



(4.9) 



(4.10) 



(4.11) 
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The only step we have still to prove is that the space where we need to 
do the orthogonalization is actually = span(l,a;, . . . ,a;"~^). This follows 
from Lemma [12] below. □ 



Lemma 12. We have the identity 



n-l 



J2 l['^et{<P^{xlx^+'))^<,j<,^, = C^ n (^i-^D- (4.12) 

xf ,l<l<k<n k=l l<i<j<n 

for some constant C„ 7^ 0. 

Proof of Lemma [TH The proof is made by induction. Assume that (14.121) is 
true for some n, which is the case for n = 2. Then we prove that (I4.12p holds 
for n + 1 . Consider 

{1 - g, x = y, 
1, y>x, (4.13) 
0, y <x. 

We will prove (I4.12p for cp instead of 0*; the statements are clearly equivalent 
by shifts 1— > + k. Set 

Then, for x > 0, G'^(x) = 1 and Gfc(x) is a polynomial of degree n — k in x 
(by evaluating the residue at z = 1), while for x < 0, G'^{x) = (because 
the residue at 00 gives zero). Therefore 

det(G'^(x;'))i<,,fe<„ = C„ n (^i-O- (4-15) 

l<i<j<n 

For the proof by induction, we apply Heine identity, 

det (G^ (xj) ) i<j,k<n det {^Pn+i (xi, yj)) i<ij<n+i 
= n! det((G^ * fn+i){yj))i<j,k<n+i- (4.16) 

The computation of the convolution leads to (G^ * {pn+i)iy) = G'k^^iy)- D 

Theorem [TT] holds for general fixed initial conditions. We want to ap- 
ply it to the alternating initial condition. For that we first have to do the 
orthogonalization with the result given in the next lemma. 
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Lemma 13. For initial conditions yj = —2j, j = 1, . . . , n, we have 

^r(^) = 7T- <P d^Tr^-T^^T, (4.17) 



and 

where, as before, p = 1 — q. In particular, $q'*(x) = 1. 

Proof of Lemma[TM The formula for is just obtained by substituting the 
initial conditions into (14.40 . Now we prove that the orthonormality relation 
(14.51) holds. For k = 0, . . . ,n — 1, the pole at w = in \Ef^ is not present, and 
for X < —2n + k, ^'^(x) = because the residue at —1 vanishes. Thus 



x) 



^ dz (b dw- 



(27ri)2 7r J^_^ zJ+i(l + pz)*-J+i (1 + w)2"-'=-i 



oo 

X 

x=—2n+k 



where we have the constraint on the integration paths |1 + -2| < |l + w|. The 
last term (the sum) equals 



-I I \ —2n+k . 
W z\ 1 + W 



1 + W / W — Z 



(4.20) 



Now the pole at if = —1 has disappeared and instead of it there is a simple 
pole sXw = z. Thus, the integral over w is just the residue dXw = z, leading 
to 

where we used the change of variable u = D 

Lemma fT3l together with Theorem [TTl leads to the kernel for the alternat- 
ing initial condition. 
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Proposition 14. For yj = —2j,j = 1,. . . ,n, the kernel K in Theorem [77] 
is given by 

+K{{ni,ti),xi;{n2,t2),X2) (4.22) 
where (p^^'^'^'^'^^'^^'^'^^^^xi, X2) is given by \2. b]) and 

Iff 

K{{ni,ti),Xi]{n2,t2),X2) = (b dw (b dz- 



(27ri)2 Jr^, (1 + 



X 



(4.23) 

i/ere Fq (resp is any simple loop, anticlockwise oriented, which includes 
the pole at z = (resp. w = —1), satisfying { — j^, z G Fq} C F_i and no 
point o/Fq lies inside F_i. 

Proof of Proposition [7^1 We substitute f l4.17p and (14.181) in the kernel (14.21) . 
Since $"'*(x) = for j < 0, we can extend the sum over k to 00. We can take 



l+pw z{l+z) 
w{l+w) 1+pz 



the sum inside the integrals if the integration paths satisfy 
Then we compute the geometric series and obtain 

(1 +z)^2+n2(l + 2;Z+pz2) 1 



< 1. 



(4.24) 

At this point both simple poles w = z and w = —{1 + z) / {1 + pz) are inside 
the integration path Fo,_i, but the integrand does not have any pole anymore 
at w = 0. Thus we will drop the in Fo,-i. Separating the contribution from 
the pole at if = z we get 



"2 



E^ni-U^i)C^U^2) = K{{n,,t,),x,;{n2,t2),X2) (4.25) 



fc=i 



27ri 



ro 



Moreover, we also have 

0*«"i'*i)'("^'*^))(xi,X2) = 0«"^'*^)'("^'*^))(xi,a;2) (4.26) 



1 f ( l+pz^ 



27ri 
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Thus the last two terms of (14.251) and fl4.26p cancel out, leading to fl4.23p . □ 

With Proposition [H] we almost obtained Theorem [H What remains to 
do is to focus far enough into the negative axis, where the influence of the 
finiteness of the number of particles is not present anymore. There the kernel 
is equal to the kernel for the initial conditions yi = —2i, i G Z. 

Proof of Theorem Ui The kernel for the fiat case is obtained by considering 
the region satisfying xi + ni + 1 < where the effect of the boundary in the 
TASEP is absent. Here the pole at w = —1 vanishes. Computing the residue 
at w = —(1 + z)/{l + pz) in Proposition [TH gives the kernel (12.51) up to a 
factor (— l)"i^"2 -^yiiich we cancel by a conjugation of the kernel. □ 



5 Proof of Theorem 3 



From Theorem [T] we have that P(x„(t) > x) = det(l — l(_oo,a;)-^l(-oo,a;))- 
We have such a situation but with x = y — sT^^^. With this change of 
variable, we get P{xn{t) > y ~ sT^^^) = det(l — 1(s^oo)-^t''''1(s,oo)) where 
^t''(6, 6) = T^^^K{xi-^iT^/^, X2-i2T^'^) (here we did not write explicitly 
the (n, t) entries). Taking into account the scaling (12.101) . we thus have to 
analyze the rescaled kernel 

i^r'(«i,ei;«2,6) (5.1) 

= T^I^K{{n{u^),t{u^)),x{u^) - (^(^2), ^(^2)), x^u^) - ^T^^') , 

with x{u) = —2n{u) + vt{u), v = 1 — y^g. In particular, we have to prove 
that, for Ui, U2 fixed, K^^'^ (or a conjugate kernel of it) converges to the ker- 
nel Ky^fC^j (kj^^mi, k^^U2, /tv^^a) uniformly on bounded sets and have 
enough control (bounds) on the decay of K'^^^'^ in the variables ^1,^2 such 
that also the Fredholm determinant converges. 

In order to have a proper limit of the kernel as T ^ 00, we have to 
consider the conjugate kernel K^^^ given by 

^ (5.2) 
The new kernel does not change the determinantal measure, being just a 
conjugation of the old one. So, in the following we will determine the limit 
of (ESI) as T ^ 00. 
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Proposition 15 (Uniform convergence on compact sets). For Ui,U2 fixed, 
according to I12.1(J\) . set 

Xi = [-2niu,) + vtiu,) - ^iT^/'] , (5.3) 
rii = n{ui), ti = t{ui). (5.4) 

Then, for any fixed L > 0, we have 

lim ir™°-'(ni,Xi;n2,X2)T^/^ = k"^/^^! (/t^ ^^^6, ^^h^^a, K'^^a) (5.5) 

T — *CXD 

uniformly for (^1,^2) ^ [~L,L]'^, with the kernel Kjs,^ given by l[2. 7| j and the 
constants k,^ and given by Ii2.13\) . 

Proof of Proposition [73 First we consider the first term in 02.41) . We thus 
consider (12.61) with the above replacements and conjugation. This term has 
to be considered only for U2 > ui. The change of variable w = —1 + 
leads then to 

Tin / ^gT2/3(go(^)-9o(^c))+Tl/3(gi(^)_gi(^,)) 

27ri /r^ z 
with Zc = {)- + ^/<^)^^ and 
go{z) = iu2 - u,)i7r'i9) + 1) ( + (1 - q)z) - (1 - v^) ln(^)) 

g^{z) = -{ul - umiT"{e) i^Hz) + ln(l - ^z)) - (6 - ^i) In(^). (5.7) 

The function go has a critical point at z = Zc- The series expansions around 

z = Zc are 

9o{z) = goizc) + {U2 - Ui)ki{z - z^f + 0{{z - z^f), (5.8) 
g^iz) = gi{z,) - (6 - 6)(1 + Vq){z - z,) + 0{{z - z,f), 



where 



vr'(^^) + l)^^ + l-7r'(^) 



(5.9) 



To prove convergence of (15. 6p we have to show that the contribution 
coming around the critical point dominates in the T — 00 limit. We do it 
by finding a steep descent pathi for go passing hy z = z^.. Consider the path 



^For an integral / — J^dze*^'^^\ we say that 7 is a steep descent path if (1) Re(/(z)) 
is maximal at some zq € 7: Re{f{z)) < Re(/(zo)) for z € 7 \ {zq} and (2) Re(/(z)) is 
monotone along 7 \ {zq} except, if 7 is closed, at a single point where Re(/) is minimal. 
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{pe*^0 G [-7r,7r)}. Then, on Tq, i^Re(ln(z)) = 0, 

^^T) n f ^^^^ \ w 7^(1 - g)psin((/.) 
-Re(ln(y^ + (1 - g).)) = (5.10) 



and 



;^He(-Ml-y.-))^-^. (5.U; 

Thus To is a steep descent path for go. Now we set p = Zc- Then, the 
real part of go{z) is maximal at z = and strictly less then g{zc) for all 
other points on Fq. Therefore, we can restrict the integration from Tq to 
Tq = {z (z Tq\\z — Zc\ < 6}. For 6 small, the error made is just of order 
0{e~'^'^^^^) with c > (c ~ 5^ as 5 ^ 1). In the integral over Fq we can use 
dii]) to get 



d^gT2/3(n2-ni)«i(2-^,)2-TV3(52_5i)(l+^)(^_2^) 



(1 + Vg)TV3 . 
2vri /rg 

We use le"^ — 1| < Ixle'"^' to control the difference between (15.121) and the 
same expression without the error terms. By taking 6 small enough and the 
change of variable {z — Zc)T^^^ = W, we obtain that this difference is just 
of order 0(T~^^^), uniformly for .^1,^2 in a bounded set. At this point we 
remain with (I5.12p without the error terms. We extend the integration path 
to Zc + ilR and this, as above, gives an error of order 0{e~^^^''^). Thus we 
have 



= 0{e~'^'^\T-^/^) (5.13) 
+ 



27ri 

Therefore, uniformly for ^1,^2 in bounded sets, 

lim (EBD = , ^ = exp i- J^^'^^i o ) (5-14) 



T^"^ y/An{u2-Ui)a \ 4(^2 -Mi)a 

with = Ki/(1 + y/qy = /t^/zth. 

Now we have to consider the second term in (12.41) . Notice that this time 
the restriction U2 > Ui does not apply. Set = — 1/(1 + ^/q). Then 

i^r(«i,6;«2,6) = ^— /^/^ ^^./3,,(.^)^^v3,,(..) (5-15) 
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with 



f,{z) = {n{9) + 9) (l-y^)ln 



1 + z 



—z 



(1 + ^ ln(l + (1 - q)z) + y/qHq) 



h{z) = {iT'{e) + 1) - ^) H-z) + V^ln((l + (1 - q)z)lq)) 

- U2{{1 - Vq) ln(l + z)- hi(l + (1 - q)z)) 
il + z)i-z) 



l + {l-q)z 



(5.16) 



f2iz) 



n"{9) 



ul{\n{l + z) + ^ln{-qz) - (1 + ^) ln(l + (1 - q)z)) 



u: 



\i^\n{l + z)+\n{-z)) 



+ 6 ln(-g2;/(l + (1 - q)z)) - ^ ln(l + z), 
f,(z) = ~H-z{l + {l-q)z)). 

The function /o has a double critical point at z 
expansions around z = z^ of the /j's are given by 

f,{z) = lK2{z-z,f + 0{iz-z,y), 

f\{z) = h{z,) + ^{U2 - u{){z - z,f + 0{{z - z,f), 

f2{z) = f2{Zc) - (6 + 6)^^(^ - ^c) + 0((^ - ^c)'), 

/3(^) = ln((l + v/g)/v/g) + 0((^ - ^.)), 
with Ki given in fl5.9p and 

(^(g)+g)(l_g)(l + ^)3 

K2 = . 



(5.17) 

Zr and the series 



(5.18) 



(5.19) 



The leading contribution in the T ^ oo limit will come from the re- 
gion around the double critical point. The first step is to choose for 70 a 
steep descent path for /q. First we consider 70 = {— pe"^,0 G [— 7r,7r)}, 
p G (0, l/(l + y^)]. The only part in Re(/o(z)) which is not constant along 70 
is the term {'7i{9)+9){l — ^) multiplied by A{z) = In |1-|--2| —a In |1-|-(1— g)z|, 
a = + v^)/(l — v^)- Simple computations lead then to 

d_ ^ sin(0)pVg(2 + - 4p(l + ^) cos(0) + p\2 - Vg)(l + ^)') 
d0 \l + z\^\l + {l-q)z\^ 

(5.20) 
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This expression is strictly less than zero along 70 except at = 0, — tt, pro- 
vided that the last term is strictly positive for 7^ 0, — tt. This is easy to 
check because the last term reaches his minimum at cos(0) = —1. Solving 
a second degree equations, we get that on p G (0, 1/(1 + y/q)) it is strictly 
positive and at p = 1/(1 + ^/q) is zero. Thus, the path 70 is steep descent 
for /o- 

But close to the critical point, the steepest descent path leaves with an an- 
gle ±7r/3. Therefore, consider for a moment 71 = {z = Zc + e~^'^^^'^^^^^^\x\, x G 
M}. By symmetry we can restrict the next computations to a; > 0. We have 
to see that B{z) = In |1 -|- z| — In |2| — a In |1 -|- (1 — q)z\ is maximum at x = 
and decreasing for a; > 0. We have 

d 

d^^^^^ ^ "2|l + z|2|^|2|l + (l-g)^|2 ^^-2^) 

(2g + 2(1 - q)^x - (1 - 3^^ + g)(l + v^) V + 2(1 - g)(l + v^) V). 



X 



The term in the second line is always positive for all x > 0. To see this, 
remark that it is a polynomial of third degree which goes to 00 as x ^ cxd 
and at X = is already positive and has positive slope. Therefore one just 
computes its stationary points and, if reals, takes the right-most one. There, 
the term under consideration turns out to be positive, which concludes the 
argument. Consequently, 71 is also a steep descent path. 

We choose a steep descent path Fq as follows. We follow 71 starting from 
the critical point until we intersect it with 70, and then we follow 79. Since Tq 
is steep descent for /o, we can integrate only on Tq = {z E Tq\\z — Zc\ < S}. 
The error made by this cut is just of order 0{e~^^) for some c = c{6) > 
(with c ~ (5^ as 5 — *• 0). Around the critical point we use the series expansions 
fl5.18p . Thus we have 

i?rK,ei;«2,6) = (^(e-^'') 

We want to cancel the error terms. The difference between (15.221) and the 
same expression without the error terms is bounded using |e^ — 1| < e'^'|x|, 
applied to x = 0{- ■ ■). Then, this error term becomes 

^ dz^-^Oiiz-z,rT,iz-z,)'T'/Mz-z,rT'^Mz-z,)) 



27ri Jr* Vq 
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for some Ci, 02,03 depending on 5. As 5 ^ 0, the Cj 1. Thus, for 6 small 
enough, we have Ci > 0. By the change of variable {z — z^jT^^^ = W we 
obtain that (15.231) is just of order 0{T^^^^). Thus we have 

i^r^K,6;«2,6) = 0(e-^^,T-i/3) (5.24) 



2vri Jrs 

The extension of the path Fq to a path going from e^^^^oo to e'^'^^^oo accounts 
into an error 0{e~'^^) only. We do the change of variable Z = K2^^T^^'^{z — Zc) 
and we define 

1/3 r- 2/3 

= (5.25) 

Then 



lim irr^(wi,6;«2,6) = / dZe3^ +("^-"^)^ '^'^ -(5i+6)^«v- 

T^oo 27ri 

(5.26) 

where 700 is any path going from e"^/^oo to e^^'^^'^oo. The proof ends by using 
the Airy function representation ( ]A.5[) . □ 

Proposition 16 (Bound for the diffusion term of the kernel). Let rii, ti, and 
Xi be defined as in Proposition [731 Then, for M2 — Mi > fixed and for any 
^1, ^2 £ the bound 



I— \ X\—X2 

(ni,ti),(n2,t2))(^3,^^3,2)2"l/3 / j ^ni-n2^-(tl-t2)/2 

1 + ) 



< const e '^^ 

(5.27) 

/ioWs /or T large enough and const independent ofT. 

Proof of Proposition UU We start with (15. 6p . The difference now is that the 
contribution coming from large |^i — ^2! can be of the same order as the one 
from go{z). We consider as path Tq = {pe^'^, (p G [— vr, tt)}. 

The difference is that now we choose p as follows. For an e with < e ^ 1 
and Zc = 1/(1 + ^), 

if 16-61 <^Ti/3, 
if6-6>^TV3, (5.28) 

By (15.101) and (15.111) . Fq is a steep descent path for go{z) plus the term 
proportional to 6 — 6 i^i fi'i(^)- So, integrating on Tq = {z = pe^'^.cj) G 



Zc + 


{6-«i)T-i/3 


(u2-mi)ki 


2;c + 


e 


(m2-'"i)ki ' 


- 


e 


(«2— ' 
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(—5, 6)} instead of Tq we do only an error of order 0{e times the value 

at = 0, for some c > 0. Thus 

LHSofdSSZD = e^''"(5°('')-s°(^^»+^'^'(5^(^)-^^(^^» (5.29) 

V ^ ^ 2m J^s z 

On Tq, the ^j-dependent term in Re{gi{z) — gi{p)) is equal to zero. With the 
same procedure as in Proposition [15] one shows that the integral is bounded 
by a constant, uniformly in T. 

It remains to estimate the first factor in (15.291) . With our choice (I5.28p . 
we need just series expansions of Qq and gi around p. Namely, by (15.81) 

T^/\9,{p) - 9,{z.)) = {u,-u,)m{p-z,fT^'\l + 0{p-z,)), 
T^"{gi{p)-gi{z,)) = (6-6)(l + v^)(p-^c)Ti/'(l + 0(p-;^c)) 

+ 0{{p-z,f)T'/\ (5.30) 

First consider the case ~ ^i| < eT^^^. We replace p given in (I5.28P into 
(I5.30p and get that the sum of the two contributions in (15.301) writes 

- ^M^^^^h + o{e) + 0{T''/')) . (5.31) 

{U2 - Ui)Ki 

0{e) comes from 0{p — z^), while the 0{T^^^^) from 0{{p — z^,)'^). Then, by 
taking e small enough and T large enough, we get 

(E3ID < -16 -61 +const. (5.32) 

In the case, 6 " 6 > £T^^^, we also replace the appropriate p given in 
(15:281) into flCTD . We explicitly use the bound eT^/^ < 6 - 6 to bound 
0((p — ZcY) < (6 — ^i)T^^^^€. Then, we obtain the following bound for the 
sum of the two contributions in (15.301) . 



16 - ^i\sT'/' f 0(T-V3) _ (1 + o{e))) < -16 - 6 



(5.33) 



by taking a fixed e small enough and then T large enough. Finally, for 
6 - 6 < eT'^^ the same result holds in a similar way. □ 

Proposition 17 (Bound for the main term of the kernel). Let rii, ti, and Xi 
he defined as in Proposition [731 Let L > fixed. Then, for given ui,U2 and 
6)6 ^ the bound 

|^^''"^(mi,6;m2,6)| < const e-(«^+«^) (5.34) 

holds for T large enough and const independent ofT. 
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Proof of Proposition \T1\ For ^1,^2 £ [~-^) is the content of Propo- 
sition [151 Thus we consider ^1,^2 £ [—L.ooY \ [— L,L]^. Define = 
(^i + 2L)T^^/'^ > 0. Then we consider a shght modification of (15. 151) . namely 

_7-l/3 y gT/o{^)+T2/3/i(2)+Tl/3/2(;2)+/3(^) 

K^^Hui,^i\U2,i2) = dz : — -7-^ (5.35) 

with /i(2) and f^iz) as in (15161) -( l5T7ll . 72(2) as /a (2) in (IHrTTll but with ^1 
and ^2 replaced by — 2L, and finally ff){z) is set to be equal to /o(-2) in (I5.16P 
plus the term 

- li ln((l + (1 - q)z)/{-qz)) - ^ ln(l + z). (5.36) 

We also chose Fq = {— pe^'^,</) G [— 7r,7r)}. In the proof of Proposition [T^ we 
already proved that Fq is a steep descent path for /o for the values p G (0, ^J- 
Also, since > 0, Re( fl5.36l) ) is also decreasing while |0| is increasing. The 
precise choice of p is 

^-^c-((el + 6)//^2)^/^ if|li + 6|<e, .33^. 
\-z,-(e/«:2)^/^ if 1^1 + 61 > 5, ^ ■ ^ 

for some small e > which can be chosen later. Let us define 

g(p) ^ gRe(T(/o(-p)-/„(2c))+T2/3(/i(-p)-/i(2,))+Tl/3(/2(_p)_/2(^^))) ^ ^g^gg^ 

Then, since Fq is a steep descent path for /o, 

_7-l/3 /• gT/o(;^)+T2/3/i(^)+rl/3/2(^)+/3(^) 



(E35D = Q{p)0{e-''^) + Q{p)^^ [ dz- 

iTXl J^s 



(5.39) 

where Fq = {—pe^'^,(f) G {—S,S)}, for a small 5 > 0. The expansion around 
= leads to 

Re(/o(-pe^'^) - /o(-p)) = -7i0'(l + 0{<j>)) (5.40) 

with 

(1 - p(l + Vg))Pv/g(l - ^){pq + (1 - p)(2 + yg)) 

2(l-p)2(l-p(l-g))2 

+ 2l(wF+(i-p(i-p))0 ^'-'^^ 

which is strictly positive for p chosen as in (15.371) . Also, Re(/i(— pe"^) — 
/i(— p)) = 720^(1 + for some bounded 72 (we do not write it down 
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explicitly since the precise formula is not relevant). Therefore, the last term 
in (I5.39P is bounded by 

const Q(p)ri/3 j' d0e-^<^"^(i+^(*))(i+^(^"''« (5.42) 

with 7 = 7i + 72T~^/^. By choosing 5 small enough and independent of T, 
and then T large enough, the error terms can be replaced by 1/2, and the 
integral is then bounded by the one on R. Thus 

(K^ < const Q(p) ^ . (5.43) 
y 7T^/3 

In the worse case, when 7-^0, which happens when p — > Zc, we have 
'j^T^/'^ ~ (^^i +^2 + ^Ly^'^ > (2L)^/^, which dominates 72 for L large enough. 
Therefore we have shown that 

|irr^(wi,ei;w2,6)| < Q{p)0{l). (5.44) 
It thus remains to find an bound on Q{p). We have, by (15.181) . 

(5.45) 

In the case |^i + ^2! < we then obtain 

Qip) < < eonste-(^^-^^) 

(5.46) 

for L ^ 1, e <^ 1. Finally, when + ^2! > we have 

Q(p) < < (5.47) 

by first choosing e > small and then T large enough. □ 

Proof of Theorem O The proof of Theorem [S] is the complete analogue of 
Theorem 2.5 in [5]. The results in Propositions 5.1,5.3,5.4, and 5.5 in [5] are 
replaced by the ones in Proposition [T51 [SI and [T71 The strategy is to write 
the Fredholm series of the expression for finite T and, by using the bounds 
in Propositions [16] and [ITJ see that it is bounded by a T-independent and 
integrable function. Once this is proven, one can exchange the sums/integrals 
and the T — > 00 limit by the theorem of dominated convergence. For details, 
see Theorem 2.5 in [5]. □ 
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6 Proof of Theorem [5 



In this section we prove Theorem [5l By Theorem [T], the right hand side of 
fl2.16p with Ui = \_{ti — Hi — Xi)/2\ can be written as Fredholm determinant 
of the kernel 

l{Xi < x,)K{{n„U),Xi; (rij^t^), Xj)l{Xj < Xj) (6.1) 

with K given in 02.41) . By the change of variable Xj = —hi + Hi + Xi, one 
obtains the Fredholm determinant of the kernel 

l{hi > Hi)K{{ni, ti), Hi + Xi- hi] {rij, tj),Hj + Xj - hj)l{hj > Hj). (6.2) 

With this preparation, we now go to the proof of Theorem 

Proof of Theorem [31 We have to analyze the kernel (12.41) with entries 

tj + Xj Hi 2tj 2Xj /r o\ 

rii = , ti = ^— , Xi = — + Hi — hi (6.3) 

and take the limit g ^ with hi, Hi fixed. The scaling of Xj might look 
different from the one in (I2.17P but, as we can see below, (16.31) with the 
last one replaced by Xj = — ^ gives the same limiting kernel. As g — 0, 
the kernel does not have a well defined limit and, as usual, we first have to 
consider a conjugate kernel. More precisely, we define 

n.i-?i2 

Jxi-X2)/2 V 



irq((xi,ti),/ii; (X2,t2),/i2) = K{{ni,ti),Xi; (ns, ^2), a;2)g^ 



q{ti-t2)/2- 



(6.4) 

What we have to prove is 

limdet(l - XHKgXH) = det(l - xhK^'^xh). (6.5) 

First we prove the pointwise convergence and then we obtain bounds allowing 
us to take the limit inside the Fredholm determinant. 

Consider the term coming from (12.61) . By the change of variable w = 
-1 + we get 

1 / d^ (v/g + (1 - g)^)*^-*^ / l-^/qz ^ "^""^ 



< . . s X , (6.6) 
2vri Jr^ z ^^i-a;2 ^^(yg + (^1 _ ' ^ ^ 

and, by inserting (16.31) . one obtains 

^ ^ |^ (v/g + (l-g)^)(l-v/g^) ^ 

7g + (1 - g)^\ ^)/v^ / + (1 - q)z\ (^-^^^/^ ^ 



{l-y/qz)z J \ {l-y/qz)z J z^^-^^^' 
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Consider q < qo for some qo < 1 fixed. Then, we can fix the path Tq 
independent of q, and the g — > hmit is easily obtained. It results in 

lim (EZD = — / ^^_e~(*^-*^)(^-^"')e(''^-''^)(^+^") (6.8) 



— Xi — ti + I2 / ^ 

where we applied f lA.4p . 

It is the turn of the term coming from fl2.5l) . We do the change of variable 
z = —w/{w + y^) and then we insert (16. 3p . The result is 

(tl+t2)/v^ / ,„ \ hi 



27ri /r„ w{\ + ^qw) \w{\ + ^qw)) \\ + ^/gu; J + 

(6.9) 



^ ^ {xi-X2)/yg+(H2-Hi)/2 



v^)(l + v^^) 

If X2 — xi > ti + t2, then for q small enough, the result is identically equal 
to zero, because the pole at = 0. If xi — X2 > ti + 12, then the result is 
also zero, because the residues at all other poles, ^/g, and 00 vanishes. 

In the other case, when |x2 — xi| < ti + 12, the apparent pole at w = — y/g 
is actually not there. So, we can choose a Fq independent of g < go for some 
go < 1. Then, we can simply take the limit g — * of the integrand, which 
leads to 

limlED = — / ^u;'^l+^2g(tl+t2)(«,-l-«,)g(x2-Xl)(«.+--l) (g_io) 



5—0 27ri /p|^ w 

= I ) 2v/(t, + t2)^-(x2-x,)^ , 

\ti + 12 — X2 + Xi y v / 

where in the last step we made the change of variable w ^ Xjw and applied 

fia. 

To have convergence of the Fredholm determinants we still need some 
bounds for large values of h\,h2. For g small enough, say g G [0,go] for 
some go < 1, we can set in (16.71) Fq = {z, \z\ = e} in the case /i2 > hi, and 
Fo = {z, \z\ = e~^} in the case /i2 < hi. Then, we get the bound 

IdnZDI < Cie-l'^^-'^^l (6.11) 

for some finite constant Ci independent of g. Moreover, in (16. 9p we can 
choose Tq = {z,\z\ = e~^}, which leads to the bound 

|(E9])| < C2e-(^2+''^) (6.12) 
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with C2 < oo independent of g < go- These two bounds are enough to have 
convergence of the Fredholm determinants. The strategy is exactly the same 
as in the proof of Theorem O □ 



7 Proof of Theorem M 

We analyze the kernel f l2.23p with the scalings 

Xi = UiT^/\ (7.1) 

t, = 7(0)T + 7' (0)«,T2/^ + ^«2tV3, (7.2) 

hi = 2U + iiT'"" (7.3) 

(See (12.241) and (I2.25P ). The strategy of he proof is the same as that for 
Theorem [3] and hence we only give the main differences. 

First we consider the first term in (12.231) . From (16.81) it is rewritten in the 
from (15.61) with gQ{z)^gi{z) replaced by 

9^{z) = {U2 - u,) (y(0)(z - l/z - 21nz) + (^ + l/z)) , (7.4) 
9i{z) = ^-^{ul - ul){z -l/z -2 Inz) - (6 - 6) In^. (7.5) 

The critical point of Qq^z) is now Zc = 1- The series expansions around z = Zc 
are 

go{z) = go{z,) + {U2 - Ui){z - z,f + 0{{z - z,f), (7.6) 
g,{z) = g,{z,) - (6 - ii){z - z,) + 0{{z - z,f). (7.7) 

The steep descent path can be taken to be Fq = {e**^, (p G [— vr, tt)}. Then the 
same arguments as in the proof of Theorem [3] give the first term in (12.71) . 
Next we consider the second term in (I2.23p . From (I6.10p it is rewritten 



in the form (I5.15P with /o(-2), /i(-z), f2{,z) replaced by 

/o(w;) = 27(0)(l/w;-w; + 21nw;), (7.8) 
fi{w) = 7'(0)(mi + U2)il/w - w + 2\nw) + {u2 - Ui){w + 1/w), (7.9) 

f2{w) = ^^{ul + ul){l/w-w + 2lnw) + {^i + ^2)lnw. (7.10) 
Their series expansions around Zc are 

fo{w) = -'^{w - If + 0{{w - If), (7.11) 

fiiw) = h{z,) + {U2 - u{){w - If + 0{{w - If), (7.12) 

f2{w) = f2{z,) + (^1 + ^2)iw - 1) + Oiiw - If). (7.13) 
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The steepest descent path is taken to be Tq = {pe**^, G [— 7r,7r)} with 
< p < 1. Using these one obtains the second term in fl2.7p . 

The bounds for the diffusion terms and the main term of the kernel are 
also proved in the same way as those of Propositions [16] and [T71 

A Some integral representations 

In this appendix we list some integral representations of the Bessel functions 
and the modified Bessel functions (we use the conventions of [1]). 

U2t) = 4-i , (A.1) 



2-Kl z 



27ri 7r„ 



27ri./ro z 



^ j Jn [2VW^) , (A.3) 



ro 



^) 4(2v^), (A.4) 



— / dt;e"'/^+"'''+^" = Ai(a2-6)exp(2aV3-a6), 



where 7oo is any path from e^'^^'^oo to e ™l'^oo. 



(A.5) 
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